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Abstract 

Based on the color-singlet model, we investigate the photoproduction of J ftp associated with a 
cc pair with all subprocesses including the direct, single-resolved, and double-resolved channels. 
The amplitude squared of these subprocesses are obtained analytically. By choosing corresponding 
parameters, we give theoretical predictions for the J/ip transverse momentum distributions both 
at the LEPII and at the future photon colliders for these subprocesses. The numerical results show 
that at the LEPII these processes can not give enough contributions to account for the experimental 
data, and it indicates that the color-octet mechanism may still be needed. At the photon collider 
with the laser back scattering photons, the resolved photon channe will dominate over the direct 
one in small and moderate pt regions with large ^/s. By measuring the J/'0 production associated 
with a cc pair, this process can be separated from the inclusive J ftp production and may provide 
a new chance to test the color-singlet contributions. 

PACS numbers: 12.38.Bx, 14.40.Lb 
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I. INTRODUCTION 



Since the discovery of J/V", heavy quarkonium has provided an ideal laboratory to inves- 
tigate the fundamental theory of strong interactions, Quantum Chromodynamics (QCD). 
Conventionally, people use the color-singlet model (CSM) [l| to describe the production 
and decay of heavy quarkonium. In order to overcome the theoretical difficulties related 
to the infrared divergences in the CSM [2, 13|] and reconcile the large discrepancy between 
the Tevatron data and the theoretical prediction given by the CSM an effective theory, 
the non-relativistic quantum chromodynamics (NRQCD) factorization formalism was pro- 
posed [sl . In NRQCD the production and decay rates of heavy quarkonium are factorized 
into the short distance parts and the long distance parts, and because the contributions of 
high Fock states are taken into account, the intermediate QQ pair that is produced in the 
short distance part can be in various states with different angular momenta and different 
colors. By introducing the color-octet mechanism (COM) in NRQCD, one may resolve the 
problem of infrared divergences in the CSM [o] and may hope to give a proper interpretation 
for the transverse momentum pt distribution of J/-ip production at the Tevatron [7|. More 
detail descriptions on many aspects of heavy quarkonium physics can be found in Ref . P 



The photoproduction of J/ip has been investigated by many authors [9|, [lOj, lUl] • In 2001 
the DELPHI Collaboration gave the measurement on inclusive photoproduction of J/ifj fl2| . 
Theoretical analysis indicates that the pt distribution predicted in the CSM is an order of 
magnitude smaller than the experimental result and the NRQCD prediction can give a good 
account for it by the COM 13|]. It has been reviewed as a strong support to the COM in 
NRQCD. The color evaporation model and the kt factorization formulism were also used to 
investigate this process Furthermore, the next-to-leading-order (NLO) QCD corrections 



to the processes 7 + 7 — ccpS'j'^^] + g and 7 + 7 ^ cc[^S'}^''] + 7 are accomplished in [l 
and the authors also give theoretical predictions for these processes at the TESLA. 

Recently, a number of studies show the importance of the heavy quark pair associated 
J/ip production in the CSM. The contributions from J/ip + c + c final states in J/'iIj inclusive 



^r3c(i)i 



production have been discussed by many authors at B factories 
the Tevatron and LHC jll], 



16 



17| and LEP 



18il, and at 



[l9|, and even been studied in the kt factorization formalism 20 1. 
Although it is a NLO process, the pt distribution can be changed and the differential cross 
section can be enhanced at large pt due to the different kinematics of the Feynman diagrams. 
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At B factories, the e"*" + e~ J/ip + c + c process gives more than half contribution to the 
total cross section of the J/ip inclusive production [2l|. In Ref. [l8|, the authors study the 
process 7 + 7 J/ip + c+c and find that the NLO process gives more contribution compared 
with that of the leading-order (LO) process 7 + 7— > J/ip + '-f a.t the LEP. In the large pt 
region, the contribution from 7 + 7— J/ip + c + c is bigger than that of the fragmentation 
process 7 + 7^0 + frag. J/ip + c + c. In Ref. j20|, the process 7 + 5'— J/ip + c + c 
was studied in the kt factorization formalism. 

In this paper, we will investigate all the subprocesses of the photoproduction of J/ip 
associated with cc in the CSM. Firstly, the full results including contributions from all 
the single and double resolved photon processes of the J/ip production associated with 
heavy quark-antiquark pair at the LEP will be presented for the first time. Secondly, these 
processes will be extended to the photon colliders and the results with different photon 
production mechanisms will be given. 

The paper is organized as follows. In section II, we give definitions of some relevant 
quantities and derive the analytical formulas of the differential cross sections for all the 
subprocesses. Numerical results are given in section III. Finally, a summary is given in 
section IV. 



II. FORMULATION AND CALCULATION 

There are three classes of subprocesses for 7 + 7— j^J/t/^ + c + c + X: As shown in 
Eq. ([1]), the direct process, where the two photons directly couple to the final heavy quarks; 
In Eq. ([2]), the single-resolved process, where one photon fluctuates to a parton (here, the 
gluon) and collide with the other photon to produce the final states; In Eq. ([3]), the double- 
resolved processes, where both the two photons fluctuate to partons to produce the final 
states. So in order to investigate the process thoroughly, the following four subprocesses 
must be calculated: 

-f + -f^J/ip + C + C (1) 

7 + ^ ^ J/^ + c + c (2) 
g + g ^ J/ip + c + c 

q + q^ J/ip + c + c (3) 
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FIG. 2: Typical Feynman diagrams for subprocesses q + q ^ J/i^ + c + c. The others can be 
obtained by reversing the fermion hnes. 



The four subprocesses involve 20, 30, 42, 7 Feynman diagrams, respectively. Fig {T] and [2] 
just show the Feynman diagrams of the processes 7 + ^? — > J/ip + c+c and q + q J/ip + c+c. 



The 
Ref. 



?eynman diagrams of the other two subprocesses are as same as those given in the 
3, [3- Following the color-singlet factorization formalism and the standard covariant 



projection method [22], the scattering amplitudes of these subprocesses can be expressed as 



Miaih) + bih) ^ ccC'+'L'^j'>){P) + cip,) + cip2)) = E E E 

Lj^Sz siS2 jk 

X {sr,S2 I SS,){LL,;SS, I JJ.){3j;3k \ 1) 
X M{a{h) + b{k2) c,(f ; si) + Cfc(f ; Ss) + c(pi) + 0(^2)), 



(4) 



where {3j;3k \ 1), {si, S2 \ SSz) and {LLz;SSz \ JJz) are the color-SU(3), spin-SU(2), and 
orbital angular momentum Clebsch-Gordan coefficients respectively for cc pairs projecting 
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out appropriate quantum numbers of the bound states. The Cl is the probabihty that 
describes a heavy quark-antiquark pair having the appropriate quantum numbers to evolved 
into a corresponding meson and, for the J/ip here, can be related to the wave function at 
the origin -R(O) or the color-singlet long distance matrix element (0|O/'''^|0) as following 

cl = im)? 

As for J/^|J production, the spin-triplet projection operator should be used, which is 
defined as 

^^^^ - ^(i;i|l^>(f;|Mf;i) 



1. 1 

2 2 



= ^ /(5,)(/^ + 2m,). 
At the same time, the color projection operator for the color-singlet state is given by 

{3j;3k\l) = 5,,/^/K. (5) 



itudes in the Feyn- 



241] to evaluate the 



We use the FeynArts |23| to generate the Feynman diagrams and amp 
man gauge, then insert the projection operators and use the FeynCalc 5 
square of the amplitudes. In calculating the subprocesses g + g J/ip + c + c, —g'^" is 
used for the polarization summation of the initial gluons and therefore the corresponding 
contribution of the ghost diagrams must be subtracted. The analytical results for every 
sbuprocesses are too tedious to be shown in this paper. In order to check the gauge invari- 
ance, the polarization vector of one initial gluon (photon) is replaced by the corresponding 
momentum in the direct and single-resolved processes and the zero results are obtained 
at the level of squared matrix element analytically. To check the gauge invariance of the 
subprocess g + g —>■ J/ip + c + c, we replace the polarization vector of one of the initial 
gluons by its momentum and use the physical polarization tensor P^^, for the polarization 
summation of the other gluon. Then the square of the amplitude vanishes. Otherwise, the 
ghost diagrams must be taken into consideration for checking the gauge invariance. Here 
the physical polarization tensor P^^ is explicitly expressed as 

P^. = -g>.u + — ^ — , (6) 

where k is the momentum of the gluon, 7] is an arbitrary light-like four vector with k -t] ^ 0. 
In the calculation, rj is set as the momentum of the other initial gluon conveniently. 



The differential cross section can be obtained by convoluting the parton level differential 
cross section with the photon density functions and the parton distribution functions of the 
photon. It is expressed as 



da{e + e — >e +e + J/ijj + c + c) = J dxidx2f-y{xi)fj{x2) 

j dxidxjfi/^{xi)f,j/^{xj)da{i + j ^ J /^j + c + c), (7) 



X 



where /^(x) is the photon density function and fi/^{x) is the parton distribution function 
of the photon. Here the labels i and j denote the parton contents of the photon, such as 
gluon and the light quarks. In the direct photon process, the distribution function f^/^{x) = 
6{l-x). 

In the photon-photon collisions, the initial photons can be generated by the 
bremsstrahlung or by the laser back scattering (LBS) from the e~^e~ collision. The spectrum 
of the bremsstrahlun g ph oton can be described by the Weizsacker- Williams approximation 
(WWA) as following ^ 



(8) 



where x = E^/Ef., a is the the fine structure constant and rrie is the electron mass. The 
definition of and Q^^^ are given by 



Qmin 1 5 (9) 

1 — X 

Qmax ~ ( 2 ^ (1 ~ 3;) + Qmini (l'^) 

where 9 is the angle between the momentum of the photon and the direction of the electron 
beam. This angle is taken as 32mrad at the LEPII. On the other hand, the laser back 
scattering can generate more energetic and luminous photons. The spectrum of the LBS 



photon is expressed as 26 1 



1 

1 — X H 4r(l — r) 

1 — X 



(11) 



where x = E^/E^, r = ^ , and the constant N is given by 



Xrrt ) ^ "''"i 2(1 -|- TTlx)'^ 
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where Xm = 4:EbEi cos^ ^. Here Eb is the energy of electron beam, Ei is the energy of the 



incident laser beam, 6 is the angle between the laser and the electron beam. The energy of 
the LBS photon is restricted by the following equation 



0<x < 



1 + X,n 



Telnov [271] argued that the optimal value of Xm is 4.83. 



(13) 




FIG. 3: The photon spectra of the WWA and LBS at ^/s = 500GeV. 



The spectra of the LBS and WWA photons are very different. While the latter depends 
only on the center-of-mass energy, the former depends on the parameter Xm also. By com- 
paring the spectra of the WWA photon at -^i = 500 GeV to the one at y/s = 1 TeV, we 
clearly see that there is no qualitative difference between them and the numerical difference 
is less than 15%. Therefore, we just show the comparison of the the spectra of WWA photon 
and that of the LBS photon at y/s = 500 GeV in the Fig. [21 And it can be seen that the 
distribution of the WWA photon is large at the small x region and tends to infinite at the 
end-point a; ~ 0. On the contrary, the distribution of the LBS photon is moderate in the 
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whole X region and get its maximum value at the largest x point. These two distributions 
can result in significant different results. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In calculating the numerical results, we choose the following parameters: Mc =1.5 GeV, 
a = 1/137, me = 0.511 MeV and the color-singlet matrix element {0\O-^^'l'{^sP)\0) = 1.4 
GeV^ 18]. The GRS99 j28| parton distribution function of photon is used and the running of 



as is evaluated by the LO formula of GRV98 29|]. Both the renormalization and factorization 
scales are fixed as ^/4:M^ + pf. The numerical results are multiplied by a factor of 1.278 to 



include the feeddown contribution from the ip' 
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FIG. 4: pt distributions of differential cross sections of J/tp+c+c production in various subprocesses 
at the LEPII. The solid, dashed, dotted and dash-dotted lines correspond to the subprocesses 
7-I-7 J/tp + c + c, ^ + g J/tp + c + c, g + g J/ip + c + c and q + q ^ J/ip + c + c, respectively. 
The experimental result of DELPHI is also presented jl^. 
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FIG. 5: pt distributions of differential cross sections of J/ip + c+c production at the photon colhder 
with the WWA photon spectrum at different ^/s. Here we use the same notations as those in FigUl 
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FIG. 6: Pt distributions of differential cross sections of J/ip + c+c production at the photon collider 
with the LBS photon spectrum at different ^/s. Here we use the same notations as those in FigHl 

The differential cross sections da /dpi for all the four subprocesses at the LEPII are shown 
in FigHl To obtain theoretical predictions, the parameters which are related to the LEPII 
experimental conditions are chosen as i/i = 197 GeV, 6max = 32 mrad and the rapidity 
cut —2 < y < 2. The constraint of center-of-mass energy for the two photons is < 35 
GeV[12']. From FigJH one can see that the direct photon subprocess is dominant at the 
LEPII with the WWA photon. The contribution from the single-resolved subprocess is 
smaller than that of the direct one by an order or more in magnitude and the contributions 
from the double-resolved subprocesses are even smaller than that of the direct one by almost 
four orders in magnitude. In the pt region that we investigated, the contribution from the 
double-resolved gluon subprocess and quark-antiquark subprocess are of the same order in 
magnitude. 
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In the future, the e~^e~ colhder may run at = SOOGeV or even at y/s = 1 TeV, and the 
LBS photon colhsion may be reahzed. Therefore, we also investigate the four subprocesses 
at ^/s = 500 GeV and ^/s = 1 TeV. For comparison, we give the theoretical results with 



the 6max of WW approximation as 20mrad [lOl] and the a™ of LBS photon as 4.83, which 
determines the maximum photon energy fraction as 0.83 [23]. In contrast to the calculation 
for the LEPII, here we do not use the constrain W < 35 GeV. 

Figj5] and Figl6] give the pt distributions of the differential cross sections at different 
center-of-mass energies with the LBS photon and WWA photon, respectively, at photon 
collider. For the WWA photon case, the direct photon production subprocess is always 
the dominant one. The contribution from the single-resolved process is less than that from 
the direct one, but larger than those from the double-resolved processes. However, in the 
case of the LBS photon, with the increase of the center-of-mass energy, the contributions 
from the single-resolved and the double-resolved gluon subprocesses are compatible with or 
even larger than that from the direct one. But the contribution from the quark-antiquark 
subprocess is much smaller than those from the other three subprocesses. ^ ^ 

FigJT] shows the parton distributions of photon in the GRS99 parametrization [28] . It 
can be seen that the gluon content is dominant in small x region and even divergent when 
X tends to zero. It is only in the large x region that the quark contents can be dominant. 

Let us first consider the subprocesses with the LBS photons as initial states. When the pt 
of J/ ip is lower or the ^/s becomes larger, the contributions from small x region partons are 
dominant. Because the LBS photon spectrum function has no singularity at small x region, 
and at the same time the gluon distribution function of photon has a great enhancement at 
small X region, the single-resolved and double-resolved gluon subprocess can be dominant 
in the J/ip production with lower pt or larger y/s. However, as for the subprocesses with 
the WWA photons as initial states, where both the photon spectrum function and the gluon 
distribution function of the photon have great enhancements at small x region, the single- 
and double-resolved subprocesses have no predominance compared with the direct one in all 
the region of pt- 

Table [H gives the integrated cross sections of every subprocesses of the photoproduction 
of J/tp associated with cc pair. From the numerical results, it can be seen that the total con- 
tributions from the resolved (including single and double resolved) subprocesses are smaller 



both the WWA photon and LBS photon at these two center-of-mass energy. Here we set 
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than that of direct one by about an order in magnitude in the case of WWA photons. On the 
contrary, in the case of LBS photons the total contributions from the resolved subprocesses 
are larger than that from the direct one for by a factor of 2.5 at ^/s = 500 GeV and 9 at 
^/s = 1000 GeV. It can also be inferred from the pt distribution presented in Fig. HI [5] and [61 
At the same time, all the integrated cross sections increase with the increase of ^/s for the 
processes initiated by the WWA photons. And in the case of LBS photons, only the cross 
section of the subprocess g + g — » J/ip + cc increased with i/i enhanced from 500GeV to 
ITeV. 

For the direct photon subprocess, our numerical result is a little different from the one in 
Ref . [isj] . The numerical results indicate that the contributions from the single-resolved and 
double-resolved processes are much less than that from the direct one at the LEPII. The 
authors of Ref. 151] have given the results of the NLO QCD corrections for the subprocesses 
7 + 7 — > J/ip + 'y and 7 + 7 ccpSg] -\-g at the TESLA. For the subprocess 7 + 7 — > J/ip + 'j, 
the K factor is smaller than one. And the QCD correction for the color-octet subprocess 
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TABLE I: Integrated cross sections of photoproduction of J /ip associated with a cc pair for different 
initial state e+e" energy ^/s and subprocesses. The pt cut of J/il) is set as pt >1 GeV. Other 
parameters and cut conditions are chosen as the one being used to calculate the pt distributions in 
the text, (units: -y/s in GeV and a in nb.) 









'^9+9 




m{wwA) 


2.06 X 10"^ 


2.91 X 10-6 


8.68 X 10-9 


6.38 X 10-9 


bm{wwA) 


3.51 X 10"^ 


2.53 X 10-5 


9.74 X lO-'^ 


2.83 X 10-s 


bm{LBS) 


5.33 X 10-"^ 


1.17 X 10^3 


2.57 X 10"^ 


1.28 X 10^5 


iom{wwA) 


4.80 X 10-^ 


5.54 X 10-5 


3.80 X 10-6 


5.22 X 10"^ 


mm{LBS) 


1.97 X 10-^ 


1.13 X 10-3 


6.66 X 10-^ 


1.21 X 10-6 



7 + 7 — > ccf^S'g] + g can enhanced the differential cross section significantly in the large pt 
region. Prom the above NLO results at the TESLA, one can expect that the NLO corrections 
to the color-singlet subprocess 7 + 7^J/'?/' + 7 could not enhance the result largely at the 
LEPII also. So the contributions from the color-octet mechanism can not be excluded in 
the inclusive J/il) photoproduction at the LEPII. The full investigation on the NLO QCD 
radiative corrections on the direct and resolved subprocesses may help us to clarify the 
situation. 

As for the photon collider with LBS initial photons, the contributions from the single- 

and double-resolved photon subprocesses become large significantly at lower and moderate 
Pt region with large \/s. This feature comes from the small x behavior of the gluon content 
distribution function of the photon, and can be checked in the future. 
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IV. SUMMARY 



In this paper, we investigate the production of J/%1) associated with a cc pair in the 
CSM in photon-photon coUisions, including the direct, single-resolved and double-resolved 
subprocesses. The formulas for the cross sections of the four subprocesses are obtained in the 
coUinear factorization formulism. Moreover, the results of the single-resolved subprocess are 
given for the first time. The numerical results show that the contributions from color-octet 
processes can not be excluded at present with the LEP experiment. 

At the photon collider with the LBS initial photons, the single-resolved and even the 
double-resolved processes will dominate over the direct one in the small and moderate pt 
regions. By measuring the final state J/ ■0 and cc pair, the process 7-I-7— > J / %l) -\- c -\- c -\- X 
can be separated from the inclusive J/ip production and could provide a channel to probe 
the parton contents of the photon. Furthermore, to separate this channel in experiment and 
compare the data with the theoretical prediction in the CSM also gives a new chance to test 
the CSM contributions. 
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